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ABSTRACT 


Subject headings: atomic processes — atomic data — line formation — X-rays: 
general 


1. Introduction 

The high spectral resolutions of the Chandra and XMM-Newton X-ray observatories 
have unveiled the useful diagnostic possibilities of oxygen K absorption. To mention a few, 
strong 0 VII and 0 VIII edges are almost ubiquitous in the spectra of Seyfert 1 galaxies 
which have been used by Lee et al. (2001) to predict of a warm dust absorber along the line 
of sight; although this conclusion has been criticized in the light of a data reanalysis (Sako et 
al. 2003), Steenbrugge et al. (2003) have detected inner-shell transitions of O III-O VI in the 
spectrum of NGC 5548 that point to a warm absorber that spans three orders of magnitude 
in ionization parameter. Moreover, Behar et al. (2003) have stressed that, in the case of both 
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Seyfert 1 and Seyfert 2 galaxies, a broad range of oxygen charge states are usually observed 
along the line of sight that must be fitted simultaneously, and may imply strong density 
gradients of 2-4 orders of magnitude over short distances. 

Schulz et al. (2002) have reported that the spectrum of the archetypical black hole 
candidate Cyg X-l shows a prominent oxygen K edge that when modelled suggests a deficient 
0 abundance approximately consistent with an abundance distribution in the interstellar 
medium (ISM). Observations of the 0 I edge and the equivalent widths of 0 II-O VIII in the 
spectrum of the black-hole binary LMC X-3 by Page et al. (2003) lead to upper limits of 
the neutral and ionized column densities that rule out a line-driven stellar wind, and thus 
suggest a black body with a multi-temperature disk. A gravitationally redshifted 0 vm Lya 
observed in absorption in the spectrum of the bursting neutron star EXO 0748-676 shows, 
unlike Fe ions, multiple components consistent with a Zeeman splitting in a magnetic field of 
around 1-2 xlO 9 G (Loeb 2003). Sako (2003) has shown that, in optically thick conditions, 
an O vui-N VII Bowen fluorescence mechanism, caused by the wavelength coincidence of the 
respective Lya (2p-ls) and Ly£ (7p-ls) doublets, can suppress O emission while enhancing 
C and N thus leading to abundance misinterpretations. 

The ls-2p resonant absorption lines from O ions with electron occupancies N < 3 are 
being used with comparable efficiency in the detection of the warm-hot intergalactic medium 
which, according to simulations of structure formation in the Universe, contains most of the 
baryons of the present epoch (Cagnoni et al. 2004). A spectral comparison of low and high 
extinction sources can lead to the separation of the ISM and instrumental components of 
the K edge (de Vries et al. 2003), the former closely resembling the edge structure of neutral 
oxygen computed by McLaughlin k Kirby (1998) with the A-matrix method (Berrington 
et al. 1987; Seaton 1987). High-resolution spectroscopy of the interstellar O K edge in X- 
ray binaries have been carried out by Juett et al. (2004), disclaiming previously proposed 
oxygen features from dust and molecular components and providing the first estimates of 
the O ionization fractions. 

' Accurate laboratory wavelength measurements for K-shell resonance lines in O VII have 
been reported by Engstrom k Litzen (1995) using laser pulses focused on compressed powder 
targets and by Beiersdorfer et al. (2003) with an electron beam ion trap; the latter method 
has been extended recently to O VI and O V (Schmidt et al. 2004). K-shell photoabsorp- 
tion of O I has been measured by Menzel et al. (1996) and (Stolte et al. 1997), and of 
O II by Kawatsura et al. (2002) where the role of theoretical data — namely level energies, 
wavelengths, (//-values, and natural widths — in experimental interpretation has been em- 
phasized. Laboratory energies for K- vacancy states in O I and O II have been obtained by 
Auger electron spectrometry (Caldwell k Krause 1993; Krause 1994; Caldwell et al. 1994), 
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but worrisome discrepancies with the photoabsorption measurements have not as yet been 
settled. 

Such experimental activity has certainly encouraged theoretical verifications which have 
employed a variety of atomic structure and close-coupling approaches. Radiative and Auger 
data related to the satellite lines of H-, He-, and Li-like oxygen have been computed in an 
approach based on the \/Z hydrogenic expansion (Vainshtein k Safranova 1971, 1978) and 
with the multiconfiguration Dirac-Fock method (Chen 1985, 1986). A theoretical study of 
the K-shell Auger spectrum of 0 I has been carried out with the multiconfiguration Hartree- 
Fock method (Saha 1994) finding good agreement with the measurements by Caldwell k 
Krause (1993). Energy positions and total and partial Auger decay rates for the [ls]2p 4 4 P 
and 2 P terms of 0 n have been determined by Petrini k Araujo (1994) with the structure 
code superstructure (Eissner et al. 1974) and photuc (Saraph 1987), the latter being 
based on an implementation of the close-coupling approximation of scattering theory. The K- 
shell photoionization cross section of the ground state of neutral oxygen has been computed 
with the .ft-matrix method by McLaughlin k Kirby (1998), giving a detailed comparison with 
the experimental energy positions for the [ls]2p 4 ( 4 P)np 3 P° and [ls]2p 4 ( 2 P)np 3 P° Rydberg 
series. More recently, the R-matrix method has also been used to compute the high-energy 
photoionization cross section of: the ground states of 0 VI in an 11-state approximation 
(Charro et al. 2000); O I taking into account core relaxation effects and the Auger smearing 
of the K edge (Gorczyca k McLaughlin 2000b); and in its relativistic Breit-Pauli mode, O I- 
O vi whereby wavelengths and /-values are obtained for the n — 2 Ka resonance complexes 
with / > 0.1 (Pradhan et al. 2003). 

We are presently interested in a realistic modelling of the oxygen K edge and in formally 
establishing its diagnostic potential. For this purpose we have computed as accurately as 
possible the required atomic data, some of which are not available for the whole O sequence, 
and used the xstar photoionization code (Kallman k Bautista 2001). Previous models 
have been fitted with two absorption edges and five Gaussian absorption lines (Juett et al. 
2004). The present data sets have been computed in a systematic approach previously used 
for the iron isonuclear sequence (Bautista et al. 2003; Palmeri et al. 2003a, b; Bautista et al. 
2004; Mendoza et al. 2004) which has allowed the determination of the efficiency of Fe K 
line emission and absorption in photoionized models (Kallman et al. 2004). 


2. Numerical methods 

The numerical approach has been fully described in Bautista et al. (2003). The atomic 

data, namely level energies, wavelengths, (//-values, radiative widths, total and partial 
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Auger widths, and total and partial photoionization cross sections are computed with a 
portfolio of publicly available atomic physics codes: autostruture, hfr, and bprm. Al- 
though relativistic corrections are expected to be negligible for the lowly ionized species of 
the O sequence, they can be more conspicuous for members with greater effective charge 
z = Z — N+l, Z being the atomic number and N the electron occupancy. For consis- 
tency, wavefunctions are calculated, unless otherwise stated, with the relativistic Breit-Pauli 
Hamiltonian 

Hbp = H nt + Hib + i?2b (1) 

where H nr is the usual non-relativistic Hamiltonian. The one-body relativistic operators 

N 

Hlb = ^2 /n( maSS ) + /n( d ) + /n(so) (2) 

n=l 

represent the spin-orbit interaction, / n ( so), the non-fine-structure mass variation, / n (mass), 
and the one-body Darwin, / n (d), correction. The two-body Breit operators are given by 

Hzb = ^ ^ 9nm{.S o) H - Qnm(s s) "h p nm (css) 4- Qnm{ b) "b 9nmi,Oo) (3) 

n<m 

where the fine-structure terms are g n m(so) (spin-other-orbit and mutual spin-orbit), p nm (ss) 
(spin-spin), and the non-fine-structure counterparts g nm (css) (spin-spin contact), g nm {d) 
(two-body Darwin), and g n m(oo ) (orbit-orbit). 

Brief outlines of the programs are given as follows specially in connection with the 
approximations used. 


2.1. AUTOSTRUCTURE 

autostructure (Badnell 1986, 1997) computes fine- structure level energies, radiative, 
and Auger rates in a Breit-Pauli framework including to order a 2 Z 4 the one- and two-body 
operators in Eq. (3), a being the fine-structure constant. This code may be regarded as an 
extension of the SUPERSTRUCTURE atomic structure package (Eissner et al. 1974) wherein 
configuration interaction (Cl) wavefunctions of the type 

*(LS) = (4) 

i 

are constructed from single-electron orbitals generated in a statistical Thomas-Fermi-Dirac 
potential and optimized variationally. Continuum wavefunctions are constructed within the 
distorted-wave approximation. Core relaxation effects (CRE) are investigated by comparing 
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an ion representation where all electron configurations have a common basis of orthogonal 
orbitals, to be referred to hereafter as approximation AS1, with one where each configuration 
has its own set, approximation AS2. 


2 . 2 . HFE 

In the hfr package (Cowan 1981), orbital bases are obtained for each configuration by 
solving for the spherically average atom the Hartree-Fock equations obtained by applying 
the variational principle to the configuration average energy. Relativistic corrections neglect 
the Breit terms listed in Eq. (3). Ab-initio eigenstates and eigenvalues are obtained and used 
to calculate A-values and wavelengths for each possible transition. Autoionization rates are 
calculated in a perturbation theory scheme where the radial functions of the initial and final 
states are optimized separately, and Cl is accounted for only in the autoionizing state. Data 
sets computed with hfr are labelled HF1. 


2 . 3 . BPRM 

The Breit-Pauli /2-matrix package (Berrington et a!. 1987; Seaton 1987) is based on the 
close-coupling approximation whereby the wavefunctions for states of an iV-electron target 
and a colliding electron with total angular momentum and parity Jn are expanded in terms 
of the target eigenfunctions 

+ T. c i*i ■ ( 5 ) 

« j 

The functions Xi are vector coupled products of the target eigenfunctions and the angu- 
lar components of the incident-electron functions; F)(r) are the radial part of the continuum 
wavefunctions that describe the motion of the scattered electron, and A is an antisymmetriza- 
tion operator. The functions $j are bound-type functions of the total system constructed 
with target orbitals. Breit-Pauli relativistic corrections have been implemented in BPRM by 
Scott & Burke (1980) and Scott & Taylor (1982), but the inclusion of the two-body terms 
(see Eq. 3) is at present in progress. Auger and radiative dampings are taken into account 
through an optical potential (Robicheaux et al. 1995; Gorczyca & Badnell 1996, 2000a) where 
the resonance energy with respect to the threshold acquires an imaginary component. 

In the present work, the A-electron targets are represented with all the fine structure 
levels within the n = 2 complex, and level energies are set to those obtained in approximation 
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AS2. Positions for states of the (N + l)-electron system are obtained from the peaks in the 
photoionization cross sections. This data set is referred to as RM 1 . 


3. Energy levels and wavelengths 

We have looked into several interactions that can influence the accuracy of the atomic 
data, namely Cl, relativistic corrections, and CRE. We have found by calculation that con- 
figuration interaction from the n = 3 complex gives rise to small contribution that can be 
generally neglected. Computations are carried out in intermediate coupling with relativistic 
corrections even though the latter are small in particular for low ionization stages. Some 
effort is focused in characterizing CRE. 

Energies calculated with approximations AS1, AS 2 , HF1, and RM1 for levels within the 
n = 2 complex of 0 ions are listed in Table 1 . Since the accuracy of the transition data (e.g. 
wavelengths and < 7 /-values) depends on the representations of both the lower (valence) and 
upper (K-vacancy) levels, detailed comparisons are carried out for each type. For the valence 
levels (see Table 2 ), a comparison of ASl and AS 2 reveals the importance of CRE which 
lower energies by as much as 20 %. Since RM 1 also neglects CRE, the best accord is with 
ASl: within 15% although the RM 1 energies are on average lower by 4%. It may be noted 
that the RM 1 calculation for 0 I has been performed in LS coupling due to the large size 
of the relativistic option. AS 2 and HFl agree to 5% and within 10% with the spectroscopic 
values (Moore 1998) which we find satisfactory. 

Regarding K-vacancy levels, the energy differences of the approximations considered 
with respect to AS2 are plotted in Fig 1 . The comparison with ASl shows that, for species 
with a half-filled L shell or with greater electron occupancies (6 < N < 8 and E < 544 eV), 
CRE again lower energies but by amounts that increase with N: A E = E(ASl) — E(AS2) < 
0.7 eV for N = 6 to A E « 2.5 eV for N = 8. On the other hand, for ions with a half-empty 
L shell or with lower occupancies (2 < N < 5 and E > 544 eV), level energies are raised 
by CRE particularly for N — 3 where A E > —1.35 eV. RM 1 behaves, as expected, in a 
similar fashion to ASl, but some inconsistencies (e.g. lowered rather than raised levels) 
are present for N < 6 perhaps due to resolution problems in determining the position of 
narrow resonances in the photoionization curves. The agreement between the AS 2 and HFl 
energies is a sound ±0.6 eV which sets a limit to the accuracy level that we can attain, and 
encourages us to adopt AS2 and HFl as our working approximations. 

AS2 and HFl K-vacancy level energies are compared with experiment and other theo- 
retical data in Table 3. It may be noted that there are no reported measurements for species 
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with 4 < TV < 6. HF1 level energies axe within 0.5 eV of the spectroscopic values for N < 3 
and the Cl values by Charro et al. (2000) for N — 2 while AS2 is undesirably higher (0.9 
eV) for the [ls]2p 3 Pf level in O VII. A convincing comparison with experiment for N > 7 
is impeded by the scatter displayed in measurements which is well outside the reported er- 
ror bars, but the photoionization measurements by Stolte et al. (1997) are 1 eV lower than 
HFl. The energy positions computed by Pradhan et al. (2003) with BPRM are noticeably 
higher than present results for N > 4, in particular for the B-like system (N = 5) where 
discrepancies as large as 5 eV are encountered. The values listed by McLaughlin & Kirby 
(1998) for O II (N = 7), which are also distinctively high (~ 3 eV), correspond to the Cl 
target they prepared for the non-relativistic A-matrix calculation of the photoionization of 
O I. They contrast with those quoted for an equivalent target by Gorczyca & McLaug hlin 
(2000b) which are in good agreement with Stolte et al. (1997). However, the values obtained 
by McLaughlin & Kirby (1998) for the [ls]2p 5 3 P° hole state in O I is in good agreement 
with AS2 and HFl. 

CRE also have an impact on the transition wavelengths where, as depicted in Fig. 2, the 
comparison of AS1 and RM1 with AS2 indicates that in general CRE lead to increasingly 
longer wavelength for N > 6 (A > 23 A) and shorter for higher charge states. It may also be 
seen in Fig. 2 that wavelengths obtained with AS2 and HFl agree to ±0.02 A setting again 
an accuracy bound. (AS2 and HFl wavelengths are also listed in Table 4.) The present 
wavelengths are compared with experiment and other theoretical predictions in Table 5. 
For N < 4 present data, in particular HFl, agree with measurements and Vainshtein & 
Safranova (1978) to ±0.02 A; in this respect, for 2 < N < 4 the values by Chen (1985, 1986) 
are somewhat long while those by Behar & Kahn (2002) are short, and for N = 3 the values 
by Pradhan et al. (2003) are long. For 5 < N < 6 the outcome is less certain due to absence 
of laboratory measurements and inconsistencies in the computed values, e.g. the order of 
the upper levels in N = 5. For N > 6 the comparison with measurements seem to indicate 
that all theoretical results are slightly short although the agreement of HFl is still within 
~ 0.02 A except for the transition from the ground level of O II and the measurement by 
Stolte et al. (1997) in O I which are just out (less than 0.05A). It must be point out again 
the undesirable scatter of the experimental data for O I. 


4. ^/-values 

In the following discussion of ^/-values, the transitions involving the strongly mixed 
[ls]2p 4 2 D 3 /2 and 2 P 3 /2 levels in O II are excluded due to extensive cancellation that result 
in completely unreliable data. Also the Is 2 J So — ► [ls]2p 3 P° intercombination transition in 
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0 VII with a very small gf -value ( gf < 10 3 ) that would require a much greater effort to 
guarantee respectful accuracy. 

A comparison of gf- values obtained with the AS1 and AS2 approximations shows that 
CRE are relatively unimportant, leading to differences much less than 10% except for transi- 
tions to the four [Is] 2s2p 2 P° relatively mixed levels of 0 vi where they jump up to ~ 15%. In 
Fig. 3 we compare the AS2 and HF1 #/- values (also listed in Table 4) showing an agreement 
within 10% except for O I where the latter are consistently 25% higher. In a comparison 
of AS2 and HF1 /-values with other theoretical data sets in Table 6, it is found that those 
by Pradhan et al. (2003) are generally 20% smaller; however, an interesting point arises 
regarding the discrepancies between AS2 and HFl where Pradhan et al. (2003) favors AS2 
for O I and HFl for the small /-value in O VI (see Table 6). 


5. Radiative and Auger widths 

The radiative width of the jth level is defined as 

Aj — Aji ( 6 ) 

i 

where A ,•* is rate (partial width) for a downward (j > i) transition. It is found by calculation 
that CRE disturb the radiative widths of K-vacancy levels in the O isonuclear sequence by 
10% or less (5% for N > 5) except when strong admixture occurs, namely for [ls]2p 4 2 D 3 / 2 
and 2 P 3/2 in 0 II and the four [ls]2s2p 2 P° levels in O VI where differences are around 20%. 

By contrast, as shown in Fig. 4, CRE reduce Auger rates by amounts that grow with N: 
from 10% for N = 4 (log A a < 14) up to 35% for N = 8; however, this trend is broken for 
N = 3 due to strong level mixing where the [ls]2s2p( 3 P°) 2 Pj are increased by a factor of 2 
(not shown in Fig. 4) even though the [ls]2s2p( 1 P°) 2 Pj levels are hardly modified (less than 
5%). The Auger widths obtained with HFl for the former levels are ~ 20% higher that AS2 
while the rest, as shown in Fig. 4, are on average lower by 10%. With this outcome we would 
estimate an accuracy rating for the Auger rates at better than 20%. As further support for 
this assertion, in Table 7 we compare our branching ratios for KLL Auger transitions with 
other theoretical and experimental values. It may be seen that the AS1, AS2, and HFl data 
are stable to better than 10%, and the accord with other data sets is for most cases better 
than 20%. 

Radiative and Auger widths computed in approximations AS2 and HFl are listed in 
Table 8. 
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6. Photoionization cross sections 

High-energy photoionizations cross sections of 0 parent ions with electron occupancies of 
3 < A r < 8 have been computed in approximation RMl; intermediate coupling has been used 
throughout except for O I because of computational size and negligible relativistic corrections. 
Following Gorczyca & Robicheaux (1999) and Gorczyca & McLaughlin (2000b), damping 
is taken into account in detail in order to bring forth the correct K threshold behavior. In 
Fig. 5 our cross sections are compared with those computed by Pradhan et al. (2003) with 
bprm in the region near the n = 2 resonances; and by Re ilman & Manson (1979) in a central 
field potential where it may be seen that, even though the background cross sections are in 
satisfactory accord, the positions and sharpness of the edges are discrepant. In the present 
approach, due to the astrophysics! interest of both the n = 2 resonances and the edge region, 
the resonance structure is treated in a unified manner showing that a sharp K edge occurs 
only in the case of O VI due to the absence of spectator Auger (KLL) channels of the type 


[ls]2p M np — » 2p M_2 np + e" (7) 

— *• [2s]2p** -1 np + e~ (8) 

— > [2s] 2 2p A *np + e - . (9) 

For the remaining ions KLL processes dominate over participator Auger (KLM) decay 

[ls]2p^np -> 2p #i_1 + e~ (10) 

-> [2s]2p M + e“ (11) 


causing edge smearing due to resonances with symmetric profiles of nearly constant width. 
Participator Auger decay, on the other hand, gives rise to the familiar Feschbach resonances — 
displayed by O vi in the region below the K threshold (see Fig. 5) — whose widths become 
narrower with n thus maintaining sharpness over the spectral head. 

In order to elucidate further the interesting properties of K resonances, we analyze the 
experimental and theoretical predictions reported for the two [ls]2p 4 ( 4 P)np and [ls]2p 4 ( 2 P)np 
Rydberg series in O I that exhibit significant discrepancies. Regarding resonance positions, 
a stringent test is to compare quantum defects, /z, defined as 

E = (12) 

where E\ im is the series-limit energy in Rydbergs; z = Z — IV + 1 is the effective charge 
(z = 1 for O i); and n is the principal quantum number. In Fig. 6, the RMl quantum 
defects for these two series are plotted relative to their respective series limits and compared 
with those obtained in photoionization measurements (Stolte et al. 1997) and the /?- matrix 



10 - 


calculation of McLaughlin & Kirby (1998). Firstly, it may be seen that the experimental 
error bars grow sharply with n; hence greater accuracy for high n is obtained by quantum- 
defect extrapolation rather than actual measurement. Secondly, the agreement of experiment 
with RM1 is very good while that with McLaughlin & Kirby (1998), who used practically 
the same numerical method (R-matrix), is only a poor 40%. This outcome may be due to 
their neglect of damping effects or their choice of a Hartree-Fock basis for the Is, 2s, and 
2p orbitals which results in a deficient representation for the K- vacancy states. Since the 
R-matrix method employs a common set of target orbitals for both valence and K-vacancy 
states, we have found the most adequate orbital basis to be that optimized variationally 
with an energy sum comprising all the states in the target model. The importance of K- 
vacancy state representations has also been stressed by Gorczyca & McLaughlin (2000b). 
Thirdly, it is found that the RM1 quantum defects (n > 3) are hardly modified by small 
target-energy shifts prior to Hamiltonian diagonalization, a standard R-matrix procedure 
to enhance accuracy by using experimental thresholds; thus, the absolute energy positions 
of these resonances mostly depends on the actual position of their series limits. Moreover, 
as shown in Fig. 7, the huge and incorrect quantum-defect decline obtained for high n in 
McLaughlin & Kirby (1998) when using two different sets of experimental thresholds, namely 
those of Krause (1994) and Stolte et al. (1997), are not due, in our opinion, to the actual 
threshold-energy shifts but to resonance misidentifications: in such cases, fi(n + 1) — / i{n ) « 
— 1 instead of approximately zero (see Fig. 7). The quantum-defect changes for n < 6 in 
McLaughlin & Kirby (1998) as thresholds are shifted are very small in accord with our 
mentioned finding. By contrast, the experimental quantum defects are sensitive to the 
series-limit positions. This can be appreciated in Fig. 7 when the quantum defects in the 
photoionization experiment of Menzel et al. (1996) are computed with the series limits by 
Krause (1994), on the one hand, and with those by Stolte et al. (1997) on the other, the 
latter resulting in unanimous agreement. 

It is inferred from the above comparison that the energy positions obtained by Auger 
electron spectrometry (Krause 1994; Caldwell et al. 1994) for the [ls]2p 4 4 P and 2 P states 
of O ii and listed in Table 3 should perhaps be scaled down by 0.4 eV. If this proposition is 
extended to the [ls]2p 5 3 P° state of O I, then the measurements by Stolte et al. (1997) and 
Krause (1994) are brought to close agreement; the only remaining experimental discrepancy 
for the latter state would then be the measurement by Menzel et al. (1996) high by 1 eV. 
Furthermore, contrary to the [ls]2p 4 ( 4 P)np 3 P° resonances of O I in our RM1 calculation, 
the absolute energy position of [ls]2p 5 3 P° is insensitive to threshold-energy shifts. This is 
a consequence of its wavefunction being mostly represented by the second term of the close- 
coupling expansion (5). Therefore, a useful parameter in comparisons with experiment — in 
as much as it gives an indication of a balanced close-coupling expansion and experimental 
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consistency — is the energy interval between the lowest n = 2 and n — 3 resonances, AE(2, 3) 
say. In Table 9, the RMl AE{2, 3) for the whole O isonuclear sequence are tabulated and 
compared with experimental and theoretical estimates. The 0.07 Ryd (1 eV) discrepancy 
in the AE(2, 3) obtained from their measurements by Stolte et al. (1997) and Menzel et al. 
(1996) in 0 I (N — 8) denotes irregular data: since the positions for the ( 4 P)3p state differ 
by only 0.07 eV the issue is then the [ls]2p 5 3 P°. The agreement of RMl with Stolte et al. 
(1997) is very good; however, if the RMl ( 4 P) threshold is improved by replacement with 
the experimental value at Halmiltonian diagonalization, then AE(2, 3) reduces to 1.00. In 
other words, the introduction of experimental thresholds improves the positions of resonances 
with n > 3 relative to the ionization threshold but not those within n — 2, as previously 
mentioned. It may be appreciated in Table 9 that this findin g is corroborated by the AE(2, 3) 
resulting from the threshold-energy shifting in McLaughlin & Kirby (1998). It is also seen 
in Table 9 that for the O isonuclear sequence AE{2, 3) as z for N > 2 and AE(2, 3) « z + 1 
for N — 1 if AE(2, 3) is expressed in Ryd. 

Regarding K-resonance widths in 0 I, from the Auger rates listed in Table 8 widths of 
159 meV (AS2) and 163 meV (HFl) are obtained for [ls]2p 5 3 P°. These are in good accord 
(~ 15%) with the measurements of 140 meV (Krause 1994; Menzel et al. 1996) and 142 
eV (Stolte et al. 1997) and the theoretical estimates of 169 meV (Saha 1994) and 139 meV 
(Petrini & Araujo 1994). The value of 185.22 meV listed by McLaughlin & Kirby (1998), on 
the other hand, is somewhat high. The latter authors also predict for the [ls]2p 4 ( 4 P)np 3 P° 
widths that decrease with n. As previously mentioned, the utter dominance of KLL Auger 
decay (see Eqs 7-9) gives rise to resonances of constant widths that have been confirmed 
with ASl for 3 < n < 5 to be at 103 meV and with a KLM branching ratio of 5 parts in 
1000 or less. 


7. Opacities 

In order to evaluate the accuracy of the present calculations, all the previously discused 

atomic data were included in the modelling code XSTAR (Kallman & Bautista 2001). The 

opacities for six different ionization parameters (0.001 < £ < 10) were obtained with a 

simple model including the following values: gas density of 10 7 * * * * 12 cm -3 , X-ray luminosity 
source of 10 44 erg/s and the solar abundances for H, He and O elements. Since the density 
remains constant, the well known (Tarter et al. 1969) form for the ionization parameter is 
implemented: 

£ = L/nR 2 


(13) 
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being L the luminosity of the source, n the density of the gas and R the distance to the 
emiting source. Thus, the monocromatic opacities are calculated assuming that ionization 
and heating are principally due to an external source of continuum photons. In Fig. (8) the 
results for each log(g) are plotted in the energy range of 400 - 800 eV, where the K edges 
for the different oxygen ions appears according to the ionization parameter. 

The top-left plot (log(£) = -3) represents the opacity for a low-ionized plasma, where a 
sharp K edge appears near to the 545 eV region, due to 0 I as well as the resonance structure 
near 530 — 540 eV that correspond to the Ka resonances of O I and O II. In the second 
graph in the left panel for log(£) = —2, the O II K edge begins to be strong in the ~ 565 
eV region, while in the bottom-left plot (log(£) = —1.5) this edge disappears when the O III 
and O IV K edges are evident at ~ 600 eV and ~ 630 eV respectively. For the top-right 
plot which correspond to log(£) = —1, the O IV contribution continues been strong, but in 
this opacity the O V K edge is predominant at ~ 665 eV. The O V edge is also the most 
important feature in the log(£) = 0 case, where the ionization of the gas is high enough to 
allow the participation of O VI in the opacity, as can be seen in the small K edge near to 700 
eV. The last plot (bottom-right) reproduces the case of a highly ionized gas with log(£) = 1, 
and the opacity is dominated by the O v, O VI and O VII K edges in the high energy region. 
In these stages, the K (3 and higher resonances apper as strong absorption features causing a 
smeared K edge in the opacities, the opposite characteristic of the low-ionized opacities. 

This behavior at the two extremes of ionization parameter is similar to that reported 
for the iron isonuclear sequence (Palmeri et al. 2002). An important difference with iron is 
that the oxygen opacities show important absoption features the region near the K edge, 
because the L-shell in the neutral oxygen ions is not completly full. As a result there are 
stronger absorption channels open in the oxygen case. 


8. Discussion and conclusions 

Extensive computations of the atomic data required for the spectral modelling of K- 
shell photoabsorption of oxygen ions have been carried out in a multi-code approach. This 
method has proven to be — as in a previous systematic treatment of the Fe isonuclear sequence 
(Bautista et al. 2003; Palmeri et al. 2003a, b; Bautista et al. 2004; Mendoza et al. 2004) — 
effective in sizing up the relevant interactions and in obtaining estimates of accuracy ratings, 
the latter often overlooked in theoretical work. Inter-complex Cl and relativistic corrections 
are found to be small, but CRE are shown to be conspicuous in most processes involving 
K- vacancy states. The inclusion of CRE in ion representations is therefore compelling which 
in practice is managed (e.g. in the hfr and AUTOSTRUCTURE codes) by considering non- 
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orthogonal orbital bases for K- vacancy and valence states; otherwise (e.g. BPRM) orthogonal 
orbital bases are optimized by minimizing energy sums that span both type of states. 

By considering several approximations and comparing with reported experimental and 
theoretical data, the present level energies and wavelength for 0 ions with electron occupan- 
cies N < 4 axe estimated to be accurate to within 0.5 eV and 0.02 A, respectively. For ions 
with N > 4, the absence of measurements in particular for 5 < IV < 6, the wide experimental 
scatter for N > 6, and discrepancies among theoretical data sets make the estimate of reli- 
able accuracy ratings less certain. It is found that for N > 4 our level energies (wavelengths) 
are probably somewhat high (short) by up to 1.5 eV (0.05 A) in neutral oxygen. Our best 
results are from approximation HF1. 

The radiative data (#/- values and total widths) computed with approximations AS2 and 
HF1 are expected to be accurate to 10% except in 0 I where the rating is not better than 
25% and in transitions involving strongly mixed levels, namely the [ls]2p 4 2 D 3 / 2 and 2 P 3 / 2 
levels in O II and to a lesser extent the four [ls]2s2p 2 Pj levels in O VI, where it is unbound. 
Admixture in the latter levels also affect Auger rates which otherwise are expected to be 
accurate to 20%. The complete and ranked set of radiative and Auger data for transitions 
involving K-vacancy levels in 0 ions is one of the main contributions of the present work. 

A second contribution are the high-energy photoabsorption cross sections for ions with 
N > 3 that depict in detail the resonance structure of the K edge. Such structure has been 
shown to be dominated by KLL damping (often neglected in previous work) which gives the 
edge a blunt imprint with arguably diagnostic potential. A comparison with the exper im ental 
K resonance positions in O I results in conclusive statements about the importance of R- 
matrix computations with well balanced close-coupling expansions, i.e. where both terms 
of expression (5) are separately convergent. In this respect, the accuracy of the energy 
interval between the n = 2 and n = 3 components of a K resonance series gives a measure 
of the close-coupling balance, and for the 0 isonuclear sequence, this energy interval when 
expressed in Ryd has been shown to be close to the effective charge, z = Z — N + 1, for 
members with N >2 and near z + 1 for N = 1. This finding could be exploited in spectral 
identification. 

As a follow-on to the present work, we intend to treat with the same systematic method- 
ology the Ne isonuclear sequence and then other sequences of astrophysical interest (Mg, Si, 
S, Ar, and Ca), always bearing in mind the scanty availability of laboratory measurements 
which as shown in the present work are unreplaceable requisites in theoretical fine-t unin g. 

MAB acknowledges partial support from FONACIT, Venezuela, under contract No. Sl- 
20011000912. TWG was supported in part by the NASA Astronomy and Physics Research 
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Fig. 1. — Energy differences for K- vacancy levels in the 0 isonuclear sequence computed 
with approximation AS2 and those with ASl (open circles), HF1 (solid circles) and RMl 
(solid triangles). Differences with ASl and RMl are due mainly to CRE. The agreement 
between AS2 and HF1 is within ±0.5 eV. 


Fig. 2. — Wavelength differences for Ka transitions in the O isoelectronic sequence computed 
with approximation AS2 and those with ASl (open circles), HFl (solid circles), and RMl 
(solid triangles). The agreement between AS2 and HFl is within ±0.02 A. 


Fig. 3. — Comparison of ^/-values for Kor transitions in O ions computed with approxima- 
tions AS2 and HFl. It may be seen that, while the agreement for most transitions is within 
10%, the HFl g /-values for O I are consistently larger by 25%. 


Fig. 4. — Comparison of Auger widths (s -1 ) for K- vacancy levels in O ions computed with 
approximations AS2 with those of ASl (filled circles) and HFl (open circles). Differences 
between ASl and AS2 (as large as 35%) are due to CRE. It may also be seen that HFl are 
on average 10% lower than AS2. 

Fig. 5. — High-energy photoionization cross sections of 0 ions showing the structure of the 
K edge. Black curve — RMl; red curve — Pradhan et al. (2003); green curve — Reilman 
k Manson (1979). 


Fig. 6. — Quantum defects for the (a) [ls]2p 4 ( 4 P)np 3 P° and (b) [ls]2p 4 ( 2 P)np 3 P° (3 < 
n < 6) resonance series of O I plotted relative to the respective threshold energy. Circles — 
present RMl results. Filled circles — experiment of Stolte et al. (1997). Filled triangles — 
R-matrix data by McLaughlin k Kirby (1998). 


Fig. 7. — Quantum defects for the (a) [ls]2p 4 ( 4 P)np 3 P° and (b) [ls]2p 4 ( 2 P)np 3 P° (3 < n < 
10) resonance series of O I plotted relative to the respective threshold energy. Circles — 
experiment of Krause (1994) and Menzel et al. (1996). Filled circles — experiment of Stolte 
et al. (1997). Triangles — R-matrix quantum defects of McLaughlin k Kirby (1998) relative 
to the experimental thresholds of Krause (1994). Filled triangles — R-matrix quant um 
defects of McLaughlin k Kirby (1998) relative to the experimental threshold by Stolte et al. 
(1997). 
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Table 1. Calculated energies (eV) for valence and K- vacancy levels in 0 ions 


N i Level AS1 AS2 HFl RM1 N i Level AS1 AS2 HFl RM1 


2 

1 

Is 2 x So 

0.000 

0.000 

0.000 


2 

2 

[ls]2p 3 P? 

566.9 

567.7 

568.1 


2 

3 

[ls]2p *PJ 

572.5 

573.5 

573.7 


3 

1 

2s “Sjy 2 

0.000 

0.000 

0.000 

0.000 

3 

2 

[ls]2s2p 2 PJ /2 

561.7 

563.0 

562.7 

563.4 

3 

3 

[ls]2s2p 2 P§ /2 

561.7 

563.1 

562.7 

563.4 

3 

4 

[ls]2s2p 2 P° /2 

567.4 

568.6 

568.0 

568.0 

3 

5 

Jls]2s2p 2 P° /2 

567.4 

568.6 

568.1 

568.0 

4 

1 

2s 2 % 

0.000 

0.000 

0.000 

0.000 

4 

2 

[ls]2p !p| 

554.2 

554.5 

554.4 

554.1 

5 

1 

2 P 2p ?/2 

0.000 

0.000 

0.000 

0.000 

5 

2 

2 P 2 P ?/2 

0.057 

0.049 

0.047 

0.047 

5 

3 

[ls]2p 2 2 D 5/2 

544.1 

544.6 

544.3 

544.4 

5 

4 

[ls]2p 2 2 D 3/2 

544.1 

544.6 

544.3 

544.4 

5 

5 

[ls]2p 2 2 P 1/2 

544.8 

545.2 

545.0 

545.3 

5 

6 

[ls]2p 2 2 P 3/2 

544.8 

545.3 

545.0 

545.3 

5 

7 

[ls]2p 2 2 S 1/2 

546.5 

547.1 

547.3 

548.2 

6 

1 

2p 2 3 P 0 

0.000 

0.000 

0.000 

0.000 

6 

2 

2p 2 3 Pj 

0.016 

0.014 

0.014 

0.014 

6 

3 

2p 2 3 p2 

0.044 

0.037 

0.040 

0.040 

6 

4 

2p 2 *D 2 

2.958 

2.750 

2.848 

2.711 

6 

5 

2p 2 x So 

5.390 

5.000 

5.237 

6.089 

6 

6 

[ls]2p 3 3 D| 

537.5 

536.9 

537.0 

536.9 

6 

7 

(ls]2p 3 3 D| 

537.5 

536.9 

537.0 

536.9 

6 

8 

[ls]2p 3 3 Df 

537.5 

536.9 

537.0 

537.0 

6 

9 

[ls]2p 3 3 Sf 

538.0 

537.4 

537.7 

538.0 

6 

10 

[ls]2p 3 3 P 2 

538.9 

538.4 

538.6 

539.4 

6 

11 

[ls]2p 3 3 P? 

538.9 

538.5 

538.7 

539.4 

6 

12 

|ls]2p 3 3 Pg 

538.9 

538.5 

538.7 

539.4 

6 

13 

[ls]2p 3 

540.8 

540.3 

540.6 

540.7 

6 

14 

[ls]2p 3 3 Pf 

542.2 

541.9 

542.2 

543.1 


7 

1 

2 P 3 4S 3/2 

0.000 

0.000 

0.000 

0.000 

7 

2 

2P 3 2d I/2 

3.972 

3.711 

3.770 

3.615 

7 

3 

2p 3 2 D| /2 

3.975 

3.714 

3.770 

3.614 

7 

4 

2p 3 2 P| /2 

5.393 

5.064 

5.131 

5.655 

7 

5 

2p 3 2 P t /a 

5.395 

5.066 

5.131 

5.654 

7 

6 

[ls]2p 4 4 P 5/2 

533.2 

531.8 

531.7 

532.5 

7 

7 

[ls]2p 4 4 P 3/2 

533.3 

531.9 

531.7 

532.5 

7 

8 

[ls]2p 4 4 P 1/2 

533.3 

531.9 

531.7 

532.5 

7 

9 

[ls]2p 4 2 D 5/2 

537.3 

536.1 

536.1 

537.0 

7 

10 

[ls]2p 4 2 D 3/2 

537.3 

536.1 

536.1 

537.0 

7 

11 

[ls]2p 4 2 P 3/2 

537.4 

536.2 

536.3 

537.5 

7 

12 

[ls]2p 4 2 P 1/2 

537.4 

536.2 

536.3 

537.5 

7 

13 

[ls]2p 4 2 S 1/2 

539.3 

538.4 

538.3 

540.4 

8 

1 

2p 4 3 P 2 

0.000 

0.000 

0.000 

0.000 

8 

2 

2p 4 3 Pi 

0.023 

0.020 

0.018 

0.000 

8 

3 

2p 4 3 P 0 

0.033 

0.029 

0.027 

0.000 

8 

4 

2p 4 1 D 2 

2.353 

2.181 

2.200 

2.149 

8 

5 

2p 4 1 S 0 

4.244 

3.975 

3.979 

4.833 

8 

6 

[ls]2p 5 3 Pg 

531.4 

528.8 

528.2 

529.7 

8 

7 

[ls]2p 5 3 Pf 

531.4 

528.8 

528.2 

529.7 

8 

8 

[ls)2p 5 3 P° 

531.4 

528.9 

528.3 

529.7 

8 

9 

[ls]2p 5 J Pf 

533.5 

531.2 

530.8 

532.7 
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Table 2. Comparison of valence-level energies (eV) 


N 

Level 

Expt a 

AS1 

AS2 

HF1 

RM1 

5 

2p 2 P? /2 

0.000 

0.000 

0.000 

0.000 

0.000 

5 

2 P 2p ?/2 

0.048 

0.057 

0.049 

0.047 

0.047 

6 

2p 2 3 P 0 

0.000 

0.000 

0.000 

0.000 

0.000 

6 

2p 2 3 Pl 

0.014 

0.016 

0.014 

0.014 

0.014 

6 

2p 2 3 P 2 

0.038 

0.044 

0.037 

0.040 

0.040 

6 

2p 2 2 D2 

2.514 

2.958 

2.750 

2.848 

2.711 

6 

2p 2 1 So 

5.354 

5.390 

5.000 

5.237 

6.089 

7 

2 P 3 4S 3/2 

0.000 

0.000 

0.000 

0.000 

0.000 

7 

2p 3 2 ^ /2 

3.324 

3.972 

3.711 

3.770 

3.615 

7 

2P 3 2 D ? /2 

3.327 

3.975 

3.714 

3.770 

3.614 

7 

2P 3 2 P ^ /2 

5.017 

5.393 

5.064 

5.131 

5.655 

7 

2p 3 2 P? /2 

5.018 

5.395 

5.066 

5.131 

5.654 

8 

2p 4 3 P 2 

0.000 

0.000 

0.000 

0.000 

0.000 

8 

2p 4 3 Pi 

0.020 

0.023 

0.020 

0.018 

0.000 

8 

2p 4 3 P 0 

0.028 

0.033 

0.029 

0.027 

0.000 

8 

2p 4 1 D 2 

1.967 

2.353 

2.181 

2.200 

2.149 

8 

2p 4 1 Sq 

4.190 

4.244 

3.975 

3.979 

4.833 


a Spectroscopic tables (Moore 1998) 


i 


Table 3. Comparison of K- vacancy level energies (eV) 


N State AS2 HFl Expt Theory 


2 

[ls]2p 3 P° 

567.7 

2 

[ls]2p *Pf 

573.5 

3 

[ls]2s2p( 3 P°) 2 P° 

563.1 

3 

[ls]2s2p( 1 P°) 2 P° 

568.6 

4 

[ls]2p ip; 

554.5 

5 

[ls]2p 2 2 D 3/2 

544.6 

5 

[ls]2p 2 2 P 

545.3 

5 

[ls]2p 2 2 S 1/2 

547.1 

6 

[ls]2p 3 3 D° 

536.9 

6 

[ls]2p 3 3 S° 

537.4 

6 

[ls]2p 3 3 P° 

538.5 

7 

[ls]2p 4 4 P 

531.9 

7 

[ls]2p 4 2 D 

536.1 

7 

[ls]2p 4 2 P 

536.2 

7 

[ls]2p 4 2 S 

538.4 

8 

[ls]2p s 3 P° 

528.8 


568.1 

568.6 a 

568. 2 e 

573.7 

574.0* 

573. 6 e 

562.7 

562.6* 

562.3 f 

568.1 

568.2* 

567.0*" 

554.4 


554.6* 

553.158 

544.3 


545.8* 

545.0 


547.0* 

547.3 


552.1* 

537.0 


537.2* 

537.7 


538.6* 

538.7 


540.7* 

531.7 

530.8(3) b 

532.9* 


530.42(4) c 

531.0* 1 

534.1* 

536.1 


535. 3 b 
538.7* 

536.3 

535.7(3) b 

535.9 h 


535.24(4) c 

539.2* 

538.3 


538. l h 
542.9* 

528.2 

527.2(3) b 

528.8* 


526.79(4) c 

528.22* 


527.85(10) d 

528.33* 


a Spectroscopic tables (Moore 1998) 

b Auger electron spectrometry (Krause 1994; Caldwell et al. 1994) 
c Photoionization experiment (Stolte et al. 1997) 
d Photoabsorption measurements (Menzel et al. 1996) 
e Non-relativistic Cl calculation (Charro et al. 2000) 
f i?-matrix calculation by Pradhan et al. (2003) 
g Multiconfiguration Dirac-Fock calculation by Chen (1985) 
h #- matrix calculation (Gorczyca & McLaughlin 2000b) 

^-matrix calculation (McLaughlin & Kirby 1998) 
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Table 4. Calculated wavelengths and ^/-values for K a transitions 


N 

t 

3 

AS2 

HF1 

N 

i 

3 

AS2 

HFl 

A 

(A) 

gfij 

A 

(A) 

gfij 

A 

(A) 

9 fij 

A 

(A) 

gUj 

2 

1 

2 

21.84 

1.22E-4 

21.82 

7.35E-5 

6 

4 

13 

23.06 

1.35E+0 

23.06 

1.35E+0 

2 

1 

3 

21.62 

8.27E-1 

21.61 

7.9GE-1 

6 

4 

14 

23.00 

5.09E-1 

22.99 

4.98E-1 

3 

1 

2 

22.02 

4.93E-1 

22.04 

4.37E-1 

6 

5 

14 

23.09 

3.37E-1 

23.09 

3.41E-1 

3 

1 

3 

22.02 

9.95E-1 

22.03 

8.81E-1 

7 

1 

6 

23.31 

4.24E-1 

23.32 

4.56E-1 

3 

1 

4 

21.80 

4.53E-2 

21.83 

4.40E-2 

7 

1 

7 

23.31 

2.83E-1 

23.32 

3.04E-1 

3 

1 

5 

21.80 

8.23E-2 

21.83 

8.09E-2 

7 

1 

8 

23.31 

1.41E-1 

23.32 

1.52E-1 

4 

1 

2 

22.36 

7.08E-1 

22.37 

6.52E-1 

7 

2 

9 

23.29 

6.22E-1 

23.29 

6.70E-1 

5 

1 

4 

22.77 

3.53E-1 

22.78 

3.33E-1 

7 

2 

10 

23.29 

1.00E-1 

23.29 

1.45E-2 

5 

1 

5 

22.74 

4.36E-1 

22.75 

4.09E-1 

7 

2 

11 

23.29 

6.72E-1 

23.28 

8.18E-1 

5 

1 

6 

22.74 

1.89E-1 

22.75 

1.80E-1 

7 

3 

9 

23.29 

4.09E-2 

23.29 

4.43E-2 

5 

1 

7 

22.66 

7.15E-2 

22.65 

6.67E-2 

7 

3 

10 

23.29 

4.76E-1 

23.29 

4.74E-1 

5 

2 

3 

22.77 

5.90E-1 

22.78 

5.61E-1 

7 

3 

11 

23.29 

1.48E-2 

23.28 

5.42E—2 

5 

2 

4 

22.77 

4.21E-2 

22.78 

4.24E-2 

7 

3 

12 

23.29 

3.98E-1 

23.28 

4.30E-1 

5 

2 

5 

22.74 

2.06E-1 

22.75 

1.95E-1 

7 

4 

9 

23.35 

2.00E-1 

23.35 

2.19E-1 

5 

2 

6 

22.74 

1.09E+0 

22.75 

1.03E+0 

7 

4 

10 

23.35 

2.04E-1 

23.35 

4.41E-2 

5 

2 

7 

22.66 

1.68E-1 

22.66 

1.52E-1 

7 

4 

11 

23.35 

1.66E-1 

23.34 

3.61E-1 

6 

1 

8 

23.09 

1.41E-1 

23.09 

1.42E-1 

7 

4 

12 

23.34 

7.48E-2 

23.34 

8.09E-2 

6 

1 

9 

23.07 

1.30E—1 

23.06 

1.30E-1 

7 

4 

13 

23.25 

2.12E-1 

23.25 

2.29E-1 

6 

1 

11 

23.03 

8.40E— 2 

23.02 

8.18E-2 

7 

5 

10 

23.35 

1.38E-2 

23.35 

9.82E-2 

6 

2 

7 

23.09 

3.17E— 1 

23.09 

3.18E-1 

7 

5 

11 

23.35 

1.67E-1 

23.34 

1.00E-1 

6 

2 

8 

23.09 

1.00E—1 

23.09 

9.95E-2 

7 

5 

12 

23.34 

1.46E-1 

23.34 

1.61E-1 

6 

2 

9 

23.07 

3.82E-1 

23.06 

3.82E-1 

7 

5 

13 

23.25 

1.01E-1 

23.25 

1.08E-1 

6 

2 

10 

23.03 

1.07E-1 

23.02 

1.05E-1 

8 

1 

6 

23.45 

4.31E-1 

23.47 

5.25E-1 

6 

2 

11 

23.03 

6.52E-2 

23.02 

6.31E-2 

8 

1 

7 

23.44 

1.44E-1 

23.47 

1.75E-1 

6 

2 

12 

23.03 

9.31E-2 

23.02 

9.12E-2 

8 

2 

6 

23.45 

1.44E-1 

23.47 

1.80E-1 

6 

3 

6 

23.10 

5.75E-1 

23.09 

5.77E-1 

8 

2 

7 

23.45 

8.61E-2 

23.47 

1.08E-1 

6 

3 

7 

23.09 

9.36E-2 

23.09 

9.44E-2 

8 

2 

8 

23.44 

1.15E-1 

23.47 

1.44E-1 

6 

3 

8 

23.09 

5.77E— 3 

23.09 

6.05E-3 

8 

3 

7 

23.45 

1.15E-1 

23.47 

1.44E-1 

6 

3 

9 

23.07 

6.11E-1 

23.06 

6.08E-1 

8 

4 

9 

23.44 

6.28E—1 

23.45 

7.85E-1 

6 

3 

10 

23.03 

3.58E-1 

23.02 

3.51E-1 

8 

5 

9 

23.52 

1.05E-1 

23.53 

1.32E-1 

6 

3 

11 

23.03 

1.31E-1 

23.02 

1.29E-1 









Table 5. Wavelength (A) comparison for Ka transitions 


N 

i 

j 

AS2 

HF1 

Expt 

Other theory 

2 

Is 2 1 Sq 

[ls]2p 3 Pf 

21.84 

21.82 

21.807* 

21.84*, 21.84> 


Is 2 1 So 

[ls]2p *P? 

21.62 

21.61 

21.6020(3) b 

21.60* 

3 

2s 2 S 

[ls]2s2p( 3 P°) 2 P° 

22.02 

22.03 

22.038®, 22.01(1)', 22.0194(16) d 

22.02', 22.06*, 22.05 k , 22.00 1 


2s 2 S 

[ls]2s2p( 1 P°) 2 P° 

21.80 

21.83 

21.82* 

21.84i,21.87 k , 21.79 1 

4 

2s 2 *S 

[ls]2p 1 P° 

22.36 

22.37 

22.374(3) d 

22.38' , 22.4V, 22.35 k , 22.33 1 
22.73 k , 22.73 1 

5 

2P 2 P ?/2 

[ls]2p 2 2 D 3/2 

22.77 

22.78 

22.74(2)' 


2P 2 P? /2 

[ls]2p 2 2 P 

22.74 

22.75 


22.67 k , 22.78 1 


2P 2 Pi/2 

[ls]2p 2 2 S 1/2 

22.66 

22.65 


22.67 k , 22. 73 1 

6 

2p 2 3 P 0 

[ls]2p 3 3 Df 

23.09 

23.09 

23.17(1)' 

23.08 k , 23. II 1 


2p 2 3 P 0 

[ls]2p 3 3 Sf 

23.07 

23.06 

23.00(1)' 

23.02 k , 23.05* 


2p 2 3 P 0 

[ls]2p 3 3 Pf 

23.03 

23.02 

" 

22.93 k , 22.98* 

7 

2p 3 4 S° 

[ls]2p 4 4 P 

23.31 

23.32 

23.36(1)® 

23.27 k , 23.30* 


2p 3 2 D° 

[ls]2p 4 2 P 

23.29 

23.28 

23.29(1)® 



2p 3 2 r° 

[ls]2p 4 2 P 

23.35 

23.34 

23.36(1)® 


8 

2p 4 3 P 

[ls]2p 5 3 P° 

23.45 

23.47 

23.51(3)®, 23.489(5) f , 23.536(2)* 
23.508(3) h 

23.45 k 


a Spectroscopic tables (Moore 1998) 
b Spectroscopic measurements (Engstrom & Litzen 1995) 
c Spectroscopy of NGC 5548 (Steenbrugge et al. 2003) 
d Electron beam ion trap measurements (Schmidt et al. 2004) 
e Auger electron spectrometry (Krause 1994; Caldwell et al. 1994) 
f Photoabsorption measurements (Menzel et al. 1996) 
s Photoionization experiment (Stolte et al. 1997) 
h lSM observations (Juett et al. 2004) 

*1/2 expansion method (Vainshtein & Safranova 1971, 1978) 
j Multiconfiguration Dirac-Fock calculation (Chen 1985, 1986) 
k R - matrix calculation (Pradhan et al. 2003) 

! HULLAC calculation (Behar &; Kalin 2002) 
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Table 6. Comparison of theoretical /-values for Ka transitions 


N 

i 

5 

AS2 

HF1 

RM2 a 

3 

2s 2 S 

(ls]2s2p( 3 P°) 2 P° 

0.744 

0.659 

0.576 

3 

2s 2 S 

[ls]2s2p( 1 P°) 2 P° 

0.077 

0.062 

0.061 

4 

2s 2 J S 

[ls]2p 1 P° 

0.708 

0.652 

0.565 

5 

2P 2 P? /2 

[ls]2p 2 2 D 3/2 

0.177 

0.167 

0.132 

5 

2P 2 P?/2 

[ls]2p 2 2 P 

0.313 

0.295 

0.252 

5 

2 P 2 P? /2 

[Is] 2p 2 2 S 1/2 

0.036 

0.033 

0.027 

6 

2p 2 3 P 0 

[ls]2p 3 3 Df 

0.141 

0.142 

0.119 

6 

2p 2 3 P 0 

[ls)2p 3 3 Sf 

0.130 

0.130 

0.102 

6 

2p 2 3 P 0 

[ls]2p 3 3 P° 

0.084 

0.082 

0.067 

7 

2p 3 4 Sf /2 

[ls]2p 4 4 P 

0.212 

0.228 

0.184 

8 

2p 4 3 P 2 

[ls]2p 5 3 P 

0.115 

0.140 

0.113 


a R - matrix data by Pradhan et al. (2003) 


Table 7. Comparison of branching ratios for KLL Auger transitions 


N 

i 

3 

ASl 

AS2 

HF1 

Expt a 

MCDF b 

CC C 

7 

[Is] 2p 4 4 P 

2p 2 

0.42 

0.42 

0.43 

0.55(13) 

0.36 

0.43 



[2s] 2p 3 

0.41 

0.40 

0.40 

0.33(3) 

0.46 

0.42 



[2s 2 ] 2p 4 

0.17 

0.18 

0.18 

0.12(2) 

0.18 

0.15 

7 

[ls]2p 4 2 P 

2P 2 

0.53 

0.53 

0.51 

0.53(9) 

0.45 

0.47 



[2s]2p 3 

0.27 

0.27 

0.28 

0.28(4) 

0.31 

0.28 



[2s 2 ] 2p 4 

0.20 

0.20 

0.21 

0.19(3) 

0.24 

0.25 

8 

[Is] 2p 5 3 P° 

2p 3 

0.53 

0.51 

0.55 

0.602 





[2s] 2p 4 

0.35 

0.35 

0.32 

0.30 





[2s 2 ] 2p 6 

0.12 

0.14 

0.13 

0.094 




a Auger electron spectrometry (Caldwell h Krause 1993', Caldwell et al. 1994) 
b MCDF calculation by M. H. Chen as quoted in Caldwell et al. (1994) 
c Close-coupling calculation (Petrini & Araujo 1994) 



Table 8. Calculated radiative and Auger widths for K-vacancy levels 


N 

i 

Level 

AS2 

HF1 

A T 

(s _1 ) 

(s- 1 ) 

A r 

(s- 1 ) 

A a 

(s- 1 ) 

2 

2 

[Is] 2p 3 P° 

6.60E+08 


3.44E+08 


2 

3 

[ls]2p J Pf 

3.93E+X2 


3.80E+12 


3 

2 

[Xs]2s2p 2 P° /2 

3.39E+X2 

6.60E+12 

3.00E+X2 

7.53E+12 

3 

3 

[ls]2s2p 2 P| /2 

3.42E+12 

5.76E+12 

3.02E+12 

6.96E+X2 

3 

4 

[ls]2s2p 2 P° /2 

3.X9E+XX 

7.56E+13 

3.12E+11 

6.67E+X3 

3 

5 

[ls]2s2p 2 P| /2 

2.90E+XX 

7.60E+X3 

2.87E+11 

6.73E+13 

4 

2 

[ls]2p !P| 

3.40E+12 

9.10E+13 

3.X0E+X2 

8.15E+13 

5 

3 

[Is]2p 2 2 D s/2 

X.34E+X2 

2.19E+14 

1.22E+12 

X.93E+X4 

5 

4 

[ls]2p 2 2 D 3 /2 

X.35E+X2 

2.18E+14 

1.28E+12 

X.92E+X4 

5 

5 

[ls]2p 2 2 P 1/2 

4.3XE+X2 

9.29E+13 

4.05E+12 

7.57E+13 

5 

6 

[ls]2p 2 2 P 3/2 

4.31E+12 

9.30E+13 

4.05E+12 

7.58E+X3 

5 

7 

[lsJ2p 2 2 S 1/2 

X.63E+X2 

1.89E+14 

1.50E+12 

1.66E+14 

6 

6 

[ls]2p 3 3 D§ 

1.05E+12 

2.37E+14 

1.06E+12 

2.XXE+X4 

6 

7 

[ls]2p 3 3 D§ 

X.05E+X2 

2.37E+X4 

1.06E+12 

2.10E+14 

6 

8 

[ls]2p 3 3 D" 

X.05E+X2 

2.37E+X4 

X.06E+X2 

2.I0E+X4 

6 

9 

[ls]2p 3 3 Sf 

4.8XE+X2 

8.87E+13 

4.8XE+X2 

7.24E+X3 

6 

10 

[ls]2p 3 3 P| 

X.20E+X2 

2.22E+X4 

X.X8E+X2 

X.98E+X4 

6 

11 

[Xs]2p 3 3 Pf 

1.20E+12 

2.22E+14 

X.X8E+X2 

X.98E+X4 

6 

12 

[ls]2p 3 3 Pg 

X.20E+X2 

2.22E+X4 

X.18E+X2 

X.98E+X4 

6 

13 

[ls]2p 3 *D§ 

3.46E+X2 

2.05E+X4 

3.46E+12 

1.84E-j-14 

6 

14 

[ls]2p 3 2 Pf 

3.62E+X2 

X.90E+X4 

3.60E+12 

1.71E+14 

7 

6 

[ls]2p 4 4 P 5/2 

8.67E+XX 

2.25E+14 

9.35E+11 

2.04E+X4 

7 

7 

[ls]2p 4 4 P 3/2 

8.67E+1X 

2.25E+14 

9.35E+XX 

2.04E+14 

7 

8 

[ls]2 P 4 4 P 1/2 

8.68E+XX 

2.25E+14 

9.35E+11 

2.04E+X4 

7 

9 

[ls]2p 4 2 D 5/2 

X.78E+X2 

2.60E+X4 

1.93E+12 

2.44E+14 

7 

10 

[ls]2p 4 2 D 3/2 

2.46E+X2 

2.36E+X4 

1.96E+12 

2.43E+14 

7 

11 

[ls]2p 4 2 P 3/2 

3.X7E+X2 

2.X0E+X4 

4.17E+12 

1.73E+14 

7 

12 

[Is] V 2 p 1/2 

3.85E+12 

1.86E+14 

4.20E+12 

X.72E+X4 

7 

13 

[Xs]2p 4 2 S 1/2 

1.96E+12 

2.43E+14 

2.X0E+X2 

2.28E+X4 

8 

6 

[ls]2p 5 3 Pg 

1.39E+12 

2.41E+14 

1.70E+12 

2.48E+X4 

8 

7 

[ls]2p 5 3 Pf 

1.39E+12 

2.41E+14 

1.70E+12 

2.48E+X4 

8 

8 

[Xs]2p 5 3 Pg 

1.39E+12 

2.41E+14 

1.70E+12 

2.48E+X4 

8 

9 

[ls]2p 5 *Pf 

2.99E+12 

2.20E+14 

3.66E+12 

2.29E+X4 
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Table 9. Comparison of AE(2,Z)* for 0 ions 


TV 

z 

P resent a 

Expt 

Other theory 

8 

1 

1.06 

1.06 b , 0.99° 

1.02® , 1.00 f 

7 

2 

1.93 



6 

3 

2.90 



5 

4 

3.73 



4 

5 

5.03 



3 

6 

5.64 



2 

7 

7.00 

7.02 d 


1 

8 

8.90 

8.90 d 



■“Energy interval in Ryd between the lowest n = 2 
and n = 3 resonances 

a Approximation RM1 for TV > 2 and AS 1 for TV < 2 

b Photoionization experiment (Stolte et al. 1997) 

c Photoabsorption measurements (Menzel et al. 
1996) 

d Spectroscopic tables (Moore 1998) 

e i?-matrix calculation (McLaughlin & Kirby 1998) 
with experimental thresholds of Krause (1994) 

f jR-matrix calculation (McLaughlin & Kirby 1998) 
with experimental thresholds of Stolte et al. (1997) 



